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In (V. Galdi et al., Phys. Rev. E57, 6470, 1998) a thorough characterization in terms of receiver 
operating characteristics (ROCs) of stochastic-resonance (SR) detectors of weak harmonic signals of 
known frequency in additive gaussian noise was given. It was shown that strobed sign-counting based 
strategies can be used to achieve a nice trade-off between performance and cost, by comparison with 
non-coherent correlators. Here we discuss the more realistic case where besides the sought signal 
(whose frequency is assumed known) further unwanted spectrally nearby signals with comparable 
amplitude are present. Rejection properties are discussed in terms of suitably defined false-alarm 
and false-dismissal probabilities for various values of interfering signal(s) strength and spectral 
separation. 

PACS numbers: 05.40.-a, 02.50.-r, 84.40.Ua 
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I. - BACKGROUND. 

Stochastic resonance (henceforth SR) is a peculiar phenomenon observed in a wide variety of physical 
systems acted by a mixture of a time-harmonic signal and (white, gaussian) noise, whereby the output 
spectral amplitude at the signal frequency shows a non-monotonic dependence i) on the noise strength at fixed 
signal frequency and amplitude, and ii) on the signal frequency, at fixed signal and noise levels. 

The possible use of SR in connection with weak signal detection experiments has been repeatedly suggested. 

The only meaningful comparison between different detection strategies is to compare the related detection 
probabilities at the same level of false alarm probabilities and available signal-to-noise ratios. In this connection, 
as demonstrated in [3| and [j], SR based detectors do not outperform matched filters. In it was further 
shown that using SR as a pre-processor (signal-to-noise ratio enhancer) does not improve the performance of a 
matched-filter detector. 

Nonetheless, strobed sign counting detectors could be interesting as computationally cheap alternatives to 
noncoherent correlators, as discussed in |3( and summarized in the next section. 



II. 



A STROBED SIGN COUNTING SR DETECTOR. 



The possibly simplest SR-paradigm is the Langevin system : 

!dV(x) t . , 
x = + Asm(uj s t + <j)) + en(t), 
ax 
x(0) = x , 

with quartic symmetric potential: 

x 2 x 4 
V(x)=-a—+b—, a,b>0, 



(2.1) 



(2.2) 



where n(t) is a stationary, zero mean, white Gaussian noise, with autocorrelation 
E[n(t)n(t + T)]=8{T). 

The probability density function (henceforth PDF) of x(t) in 12.1fl . denoted as p(x,t) is ruled by the 
Fokker-Planck equation: 



dp{x,t) d 
dt dx 



dV(x) 



dx 



A sin(w s t + ' 



p{x,t) \ + 



e 2 d 2 p{x,t) 
2 dx 2 ' 



(2.3) 



t p(a;,0) = 8(x-x ). 
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The solution of (|2.3|) . in the absence of a signal (A — 0) is an even function of x 0j- In the presence of a signal, 
even-symmetry is broken, and the asymmetry is maximum at 

f = i* -^T 1 f — ^— 7T- , fc = 0,l,2,...,2iV-l (2.4) 

where ^ is a (known Q) phase-lag introduced by the SR processor. 

Symmetry breaking of the output PDF is perhaps Q the most natural signature of A ^= in (|2.1(l . In Q 
we gave a thorough evaluation in terms of receiver (detection) operating characteristics (henceforth ROCs) of 
the possibly simplest non-parametric detector based on the above symmetry-breaking, where from the output 
samples (|2.4|) . one forms the time series 

x k = {-) k x(t k ), (2.5) 

and compares 

2N-1 

N+=J2 U ( x *) ( 2 - 6 ) 

k=0 

where [/(■) is Heaviside's step function, to a suitable threshold T. The above will be henceforth referred to as 
strobed sign-counting stochastic resonance detector (SSC-SRD). 

The SSC-SRD performance is described by the following false-alarm and false dismissal probabilities [3j: 

a = Prob{iV + > r|no signal} = I 1/2 (T + 1, 2N - T), (2.7) 



= Prob{AT+ < r|signal} = 1 - I P+ (T + 1, 2N — T), (2.8) 

where P + = Prob(a:fc > 0) and I p (x,y) is the incomplete Beta function . The related ROCs are typically worse 
by w 3dB as compared to those of the matched filter Q . 

Generalization to the case where the initial phase of the sought signal is unknown is straightforward, by 
letting: 



+ (2.9) 

where T s = 2%/u) s . The resulting unknown-initial-phase detector for A^ s ~ 10 has nearly the same performance 
as (|2.fi|l . which applies to the coherent (known initial phase) case 0, and is accordingly comparable to that of 
the noncoherent correlator (std. optimum benchmark detector for signals with unknown initial phase). 

On the other hand, the SSC-SRD is computationally extremely cheap, requiring only binary and/or integer 
arithmetics, and thus quite appealing. 

The obvious question is now related to the rejection properties of the above detector, namely to its ability to 
discriminate between spectrally nearby signals (sought and unwanted) with comparable amplitudes. 

Note in this connection that the frequency response of the SSC-SRD depends very little on the stochastic reso- 
nance condition, as shown in Fig. 1 (see §), where the steady-state value of maxt [-P+ (i)] = max t [Prob(a;(i) > 0)] 
is displayed for several values of the SNR as a function of ui s = uj s Tk, where u s is the signal angular frequency 
and Tj, = v / 27rexp(2Vo)/a is the Kramers time. In Fig. 1 Vq = a 2 /{Abe 2 ) is the normalized potential-barrier 
height. 

A numerical investigation of this issue will be the subject of the next Section. 



N+ = max U 
me(0,JV,) f-^ 



III. - ROCS IN THE PRESENCE OF NEARBY SIGNALS. 



In order to evaluate the performance of the above described SSC-SRD in the presence of spectrally nearby 
signals with comparable amplitude, we introduce the following dimensionless parameters: 

A u 

Ps = -j~> 1 = 



Acj) 



(3.1) 
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representing the unwanted-to-wanted (sought) signal-to-signal ratio (SSR), the (scaled) frequency difference, 
and the phase lag, where the suffixes w and u refer to the wanted and unwanted signal, respectively. More or 
less obviously, we also define the false-alarm (a) and false-dismissal (f3) probabilities as follows: 

a = Prob{iV + > T\s w (t) = and s u (t) ^ 0} , (3.2) 

p = Prob{iV + < T\s w (t) and s u (t) ^ 0} . (3.3) 

Representative numerical simulations are accordingly summarized in Fig.s 2 to 5 below. 

Figure 2 describes a situation where SNR = 3, and the unwanted signal frequency is pretty close to the 
sought one (7 = 1CP 4 ). Only two detection characteristics are displayed, corresponding to ps = 2, 3 for the sake 
of readability, together with the limiting curve corresponding to the absence of the unwanted signal (dashed 
line). It is seen that the detector's performance is only slightly deteriorated due to the presence of the unwanted 
signal. 

In Fig. 3 the SSR is fixed at (ps — 3), and two different (normalised) frequency separations 7 = 1CU 2 , 10~ 4 
are considered. The limiting curve corresponding to the absence of the unwanted signal (dashed line) is also 
included. Again, the detector's performance is not appreciably spoiled. 

In Fig. 4 the unwanted signal amplitude and spectral separation are fixed, ps = 3 and 7 = 1CP 4 , and two 
values of the phase lag are considered, A</> = 0, ir. Once more, the detector's performance is negligibly affected. 

Finally, Fig. 5 shows the key role played by the number of strobed samples used in i|2.9[) at fixed ps = 3, 
7 = 10~ 4 , Acf) = on the SSC-SR detector's rejection performance. It is seen that the unwanted signal is 
rejected for N > 10 3 , for which the time series length corresponds to the spectral width required to separate the 
wanted and unwanted signals. This shows that the rejection properties are essentially related to the spectral 
filtering inherent to the strobing process. 

IV. - CONCLUSIONS AND RECOMMENDATIONS. 

Numerical simulations suggest that strobed-sign-counting stochastic-resonance-based detectors besides 
being computationally cheap, and nearly as much performing as the standard non-coherent correlator, do 
display nice rejection properties in the presence of spectrally nearby signals with comparable amplitudes. 

Further study is under way to evaluate the potential of SSC-SR detectors in connection, e.g., with the search 
of weak nearby quasi-monochromatic signals in the context of the search of gravitational waves. 
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Captions to the Figures 

Fig. 1 - Frequency response of SSC-SR detector for several values of SNR. Dashed curves refer to the adiabatic 
approximation (quoted from |3j). 



Fig. 2 - ROCs of SSC-SR detector. SNR = 3,7 = 1(T 4 , Ac/) = 0, N = 100. The p s = curve is shown dashed. 



Fig. 3 - ROCs of SSC-SR detector. SNR = 3, p s = 3, A0 = 0, 7 = 10~ 2 , 10~ 4 , N = 100. The p s = curve is 
shown dashed. 



Fig. 4 - ROCs of SSC-SR detector. SNR = 3, p s = 3, 7 = 10~ 4 , = 0,tt, N = 100. The p s = curve is 
shown dashed. 



Fig. 5 - ROCs of SSC-SR detector. SNR = 3, p s = 3, 7 = 10~ 4 , A<f> = 0, N = 50, 100, 1000. The p s = curve 
is shown dashed. Close-up in the inset. 
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